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Abstract 

Every relativistic particle has 4-speed equal to c, since gr^ji/^^^r = (? ■ With the choice of 
k = +1 in the FRW metric, the cosmological scale factor a{t) has the natural interpretation of the 
radius of the sphere S'a = {x : {x,x) = a?}. Thus, a particle at rest in the cosmological frame 

has 4-speed equal to This leads us to infer that a — c, which respresents a simple kinematic 
constraint linking the speed of light to the cosmological scale factor. This drastically changes the 
fe = +1 picture from a closed deaccelerating universe to an open accelerating universe, settles the 
horizon problem, and provides for a new cosmological model more appealing to our natural intuition. 
In this paper we shall consider ramifications of this model. 
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Part I 

Discussion and Overture 

In this part of the paper, we will introduce our hypothesis, which will give us a simple kinematic constraint 
between the evolution of the scale factor and the speed of light. We will then consider various elements 
of cosmological theory that need modified if one is to consider varluminopic cosmological models. This 
will set the stage for our discussion of the dynamics in the upcoming parts of the paper. 



1 Introduction 

In subsection 1, we will review some empirical evidence which points to the k ^ +1 FRW Model. In 
subsection 2, wc will introduce a kinematic constraint on this model which will change its resulting 
dynamics. In subsection 3, wc will give an outline for the organization of this paper. 

1.1 The Appeal of the k = +1 Cosmologies 
1.1.1 Expansion 

Cosmological observation shows us that our universe is expanding. This was first discovered by Edwin 
Hubble in 1929. As geometry codifies the interplay between gravitation and matter, this is a natural 
implication of general relativity. It can be mathematically described independent of your choice of 
universe, by including a time-varying scale factor in the space-time metric. However, in the k = +1 
picture, expansion takes a very natural and elegent form. The universe takes the form of a three- 
dimensional expanding spherical surface in a four dimensional space, the radial component being time- 
like. In this picture, expansion is no more surprising than observing several dots on a balloon mutually 
moving farther apart as the balloon expands. 



2 



1.1.2 The Horizon Problem 



Cosmologist have measured to temperature of the CMBR to be 2.728 K with anisotropics in temperature 
on the level of 10~^. The problem is that this is uniform in every direction in the sky. In particular, if we 
look at two points separated by 180°, we record the same background temperature. The radiation from 
each of these points has travelled approximately 98% of the horizon distance, from the last scattering 
surface, to reach us. Therefore, the two points where the radiation originitated are seperated by a 
distance of 196% of the horizon distance. If the universe is taken to be flat, we arrive at a fundamental 
paradox. How could two separate points in space-time, not causally connected, be in thermal equilibrium? 
However, if the universe is spherical, the paradox is resolved. All points in space are causally connected, 
as the entire universe originates from a single point in the past. Thus, by virtue of the geometry, our 
model solves the Horizon Problem without the need of inflation. 



1.1.3 The Simultaneity Problem 

Another cosmological problem is that the standard Big Bang model, when set in a fc = flat FRW 
universe, violates the principle of simultaneity. In this model, the Big Bang occurs everywhere simul- 
taneously. That is, on a space time diagram, there is a horizontal cutoff somewhere. Thus for every 
observer in motion relative to the cosmolgical frame (i.e. us), there are regions of space (on the surface 
t = now, where t is coordinate time) where the big bang is currently going on, and there are regions of 
space where the universe has yet to be born. This problem is irrelevant in the spherical k = +1 FRW 
universe. 



1.2 A New Kinematic Constraint 



Despite these observations, the k = +1 FRW universe has been found unsatisfactory and inadequate 
by the cosmology community over the last several decades. We would like to introduce a fundamental 
constraint on the motion which radically changes the resulting dynamics of the fc = +1 FRW picture. 
Taking the metric to have signature (+1, —1, —1, —1), it is a basic result of relativity theory that 



dx^" dx" 
dr dr 



(1) 



For a particle at rest in the cosmological frame, and choosing a comoving coordinate system, the particle's 
four velocity must necessarily take the form 



dxi" 



dr 



(c, 0,0,0) 



(2) 



in order for the condition (QJ to be satisfied. 

Now consider the fc = +1 FRW picture. In spherical coordinates, the particle's position is (a, p, 6*, (p). 
If the particle is at rest in the cosmological frame, then p = = = 0, and the particle's four-velocity 
becomes 



dx^' 



dr 



da \ / da 

-,0,0,0) = (-,0,0,0 



(3) 



where t is the proper age of the universe (i.e. cosmologicaly time). 

Equating (O with JSJ, we have the following fundamental cosmological hypothesis: 

Hypothesis 1 The universe is spherical and the speed of light varies with cosmological time t subject to 
the constraint 

ait) = cit) (4) 

Having a speed of light which varies with cosmological time is no stranger than having a scale factor 
which varies with cosmological time. 
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1.3 Outline 



In Part HJ wc will devote our attention to necessary modifications of classical cosmology that should 
be taken into account when dealing with varluminopic theories. In we will state the metric Ansatz 
connected with our hypothesis Q). In SjSl wc will make a brief comment on physical constants. In 21 
we will show that an extra term appears in the geodesic equations when dealing with a varluminopic 
theory. This arises as the geodesic equations are geometric and not dynamical equations, and thus must 
be written with derivatives with respect to arclength ds. In converting to derivatives with respect to 
proper time, the chain rule now takes effect, thus modifying the standard geodesic equations. 

In Partini we will apply Einstein's Field Equations to our metric Ansatz, and look at the resulting 
dynamics. Thus in Sj^l we form all of the tensor structures related to the metric, such as the Ricci Tensor, 
Einstein Tensor, Christoffel Symbols, and Ricci Scalar. In tl6.2l we write the resulting field equations, 
which, unlike their isoluminopic FRW counterparts, do not yield an acceleration equation, rather pressure 
and density equations. An additional constraint is needed to specify the resulting dynamics. However, 
it uniquely establishes the equation of state in cosmological time, which was not previously the case. In 
^we find the resulting covariant divergence of the field equations, which contain additional terms due 
to the varying speed of light. One now has the freedom to constrain the total physical energy of the 
universe to be conserved, which was not previously possible. We will explore this scenerio, the so-called 
isoergic case, in Alternatively, we can impose an adiabatic condition, tantamount to constraining 
the covariant divergence of the stress-energy tensor to vanish. In this case, one is forced to pick up an 
additional dynamic variable, either G or A or both. Thus one would require an additional constraint to 
determine the resulting dynamics. 

In Part IIIII we introduce a new set of Varluminopic Field Equations. Wc start from the hypothesis 
that the varluminopic and isoluminopic gravitational actions coincide. We then derive the resulting field 
equations, in parallel to their classical derivation, while including varluminopic effects. In Jj^l wc show 
that the standard Lagrangian density can be recast in a form dependcnd only on the metric and its 
first derivative, as is done classically, even in the varluminopic case. In ^101 we derive the resulting 
field equations. Finally, in we impose the additional energy conservation constraint V^Tj^ = 0, and 
determine the resulting dynamical evolution of the universe. We note that the extra terms which appear 
in the Varluminopic Field Equations, as they appear in HIOI save us from having to add additional 
dynamical variables and constraints on the motion, as is the case when one attempts to impose this 
constraint using Einstein's Classical Field Equations. 

We would additionaly like to state that some work has been done in varluminopic cosmoliges (com- 
monly known as VSL or variablc-spced-of-light theories). See Basset et al |5j and the references contained 
therein. Barrow ^ and Magueijo 0, for instance. This paper will present a different approach on var- 
luminopic theories. 



2 Change of Variables 

The FRW metric with fc = +1 is equivalent to the form: 

ds^ ^ c^dT^ = c^dt^ - a{t)^{dp^ + sm^(p)de^ + sm^(p) sm^{e)d(f^} (5) 

the spatial component of which represents the restriction of the standard Euclidean metric in to 
the three-sphere S^. Our paper diverges from the classical analysis of this problem by imposing the 
kinematic constraint 

da , , 
M = 

We now introduce the change of variables t — > a{t) defined by: 

a{t) = I c(s) ds (7) 

Since c{t) > 0, a{t) represents a monotone, increasing function of the cosmological (coordinate) time t. 
This transformation can therefore be inverted to express: 

c(a) := c{t{a)) (8) 



4 



Under this change of coordinates, our metric Ansatz becomes: 

ds^ = c{afdT^ = da^ - a^{dp^ + sm^{p)de^ + sin2(p) sm^ {e)d4)'^} 



(9) 



3 A Note on Fundamental Physical Constants 

In Einstein's derivation of the classical field equations, he supposes that the Einstein tensor is proportional 
to the stress-energy tensor of the field, i.e. 

G^u = aT^, (10) 

where a works out to be 

SttG 

(11) 

where G is Newton's constant, which appears in Newton's Law of Gravitation: 

GMm , 

(12) 

When considering varluminopic theories, we must discover whether it is cr or G which is to be held 
constant. Had Einstein come up with the field equations 1)1 n|l first, then Newton's Law would have been 
written 

^ ac^Mm , „, 

and it would have been a which would have been considered the fundamental constant of gravitation. 

In this paper, we will take G to be the fundamental constant. We would like to remark that this is 
not entirely obvious, and is an assumption. There is no distinction in the isoluminopic models, however 
it now makes a decided difference when wc treat the speed of light as a dynamic variable. 



4 Modified Geodesic Equations 

If one considers the motion of a particle at rest in the cosmological frame (where ds^ ~ da?'), one would 
expect that the dynamics of the scale factor should be determined by solving the geodesic equations 
associated with the metric @. One might naively mistake this evolution to be nontrival, i.e. a = 0. 

.2 

This conclusion would be invalid. What is actually being observed is the obvious condition that = 0. 

In our theory; in fact, in any variable speed of light (VSL) theory, one can modify the geodesic 
equations to give dynamical equations of motion, as opposed to the standard geometric equations of 
motion. This distinction does not exist in classical general relativity, where the speed of light is a 
constant and can be moved freely through differential operators. 

The geometric geodesic equations, which are valid even in the VSL case, are the standard geodesic 
equations of motion: 

d^x^' „ dx'' dx^ , , 

ds^ ds ds ^ ' 

However, using the relation 

ds^ = citfdT"^ (15) 



where t is the local cosmological time, one sees that: 

dxi" _ dxf" 1 
ds dr c{t) 

d^a;^ _ d^T^ 1 fixi^ 1 dc{t) dt dr 

ds^ ^ dr2 c(t)2 ^ dr (c(t))2 dt d^ds 
_ d^T^ 1 dx^' 1 dc{t) dt 

dr^ c(<)2 ^ dr (c(i))3 dt d^ 



(16) 
(17) 
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Substituting into the geometric geodesic equation (|14|l . one obtains: 



rf^x^ dx" dx^ _ dx^" c'{t) dt 

dr^ ^ "^Ihlh ~ dT c{t) 'ch 



Or, alternatively, 



d^o:'' ^ dx" dx^ dxi" dc(a) dt , , 

dr2 dr dr dr da dT ^ ' 

We note that if a particle is at rest in the cosmological frame, then dt — dr. Substituting the Christoffel 
Symbols from into the modified geodesic equations yield no information on the dynamics of the scale 
factor, i.e. they are trivially satisfied. To determine the dynamics of the scale factor, we now turn to 
Einstein's Field Equations. 



5 Redshifts and Distances 

We derive the appropriate modifications for redshifts and luminosity distances, taking into account the 
changing speed of light. See Bergstrom [3] and Breton among many others, for a review of the 
isoluminopic derivation. 

5.1 Redshifts 

Light waves travel along nuU-geodcsics. Setting the 6 and cp angular coordinate displacements to zero, 
we have, for photons travelling along constant p world lines, the following: 

a^dt^ - a^dp^ = (20) 

For a photon emitted at tg and absorbed at to, we thus have 

a[t) Jo 

Suppose that the next wave peak is emitted and aborbed Ste and Sto later, respectively. Since they both 
travel through an angular displacement p, we have 

' ^ ' dt= dt (22) 

From which it becomes clear that 



ait) 7te+5te ait) 



ait) , /■*°+*'*° d(i) , 

^dt= ^dt 23 

ait) Jt^ ait) 

Assuming that d{t) and a{t) are approximately constant over the durations Stg and 5to, we find that: 

6to _ d(<e)Ao _ dite) aito) 



Ste d(to)Ae aito) aite) 
where A is the wavelength. Therefore the cosmological redshift z is given by 



(24) 



Ae a[te) 



which is the regular formula. 
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5.2 The Luminosity Distance 



The so-called luminosity distance is measured indirectly by measuring the arrival power flux T of 
light from distant objects with known intrinsic luminosities £, by the relation 



C 



(26) 



The luminosity distance is the distance we would be from the object if the universe was flat and ex- 
pansionless. To measure the expansion of the universe, it is standard practice to derive a theoretical 
relationship between the luminosity distance and redshifts of incoming light from distant sources, and 
then compare this relationship to observation. As usual, we must tweek this theoretical relationship to 
find the appropriate one for our new cosmology. 

Suppose a photon with energy is emitted from a distant source at time and with wavelength 
Ag. Thus 

ha{te) 



Its observed energy at arrival will be 

ha{to) a{to) \e 



d{to) 



1 



-E, 



Xo a{te) Xo a{te) (1 + z) 

The object's intrinsic luminosity C is the power of a burst of photons emitted at time te'- 

SE, 



C = 



Recalling H24(l . we see that the object's power at arrival is given by 

2 



Co 



SEq 



d{to) 



a{te) 



1 



{1 + zy 



(27) 



(28) 



(29) 



(30) 



Suppose the constant angular distance to the source is r (where we've replaced p in to avoid confusion 
with the energy density). Then the arrival flux is given by 



C 



47ra(to)2r2(l + z) 



d(to) 



d{te 



(31) 



As light travels along a null geodesic from when it is emitted at the source to to when it arrives at the 
detectors, we have 

= a{t) dt - a{t)dr (32) 

Hence 

' a{to) 



df 



= ait) 



In 



a{te) 



(33) 



This is the simple relation between the angular displacement of a distant object {a{to)r being its current 
instantaneous distance) and the scale factor. Comparing (|26|l to l|31|l . we have for the luminosity distance: 



a{to)'^ d{te) , 
«L = — 7-^ . ,. m 



a{te) 



(34) 



a(ie) d{to) 

(|25|l combined with H34|l provide a theoretical relation between the redshift and luminosity distance of 
detectable light. If one can find explicit solutions for a{t) and a(t), one can then determine the theoretical 
relationship between d^ and z and compare with known data. 

In practice, one also computes the distance modulus, defined by 



= 5 logio 

and plots this measure of distance vs. redshift. 



1 Mpc J 



25 



(35) 
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Part II 

The Classical Field Approach 

In this part, we insert our metric Ansatz ^ into Einstein's Field Equations and determine the resulting 
relationships. For the standard FRW metric, the Field Equations produce an equation for the density 
(the Fricdmann Equation) and an acceleration equation. Classically, one must in addition impose an 
equation of state to determine the dynamical outcome. One can gain insight as to how the universe 
should evolve with various equations of states, but the equation of state as a function of cosmological 
time has not been determined. 

On the other hand, we will show that our metric Ansatz instead produces a density equation and 
a pressure equation, so that one can completely determine the equation of state as a function of the 
dynamic variables a{t) and c(i). In order to get the dynamics, one must pose an additional constraint 
on the system. At this point we will distinguish two possibilities. One may take a naive approach and 
constrain the total physical energy of the universe to be constant. In this case, our so-called "isoergic" 
model, one can determine the resulting dynamics of the system. This approach cannot be done in the 
isoluminopic models (as the covariant divergence of the Field Equations vanish), so it has not been 
seriously considered. However, its philosophical implication is that, for it to be correct, there can be no 
transfer of energy from the matter content to the geometry (i.e. gravitational field). The second, more 
standard approach, is to set the covariant divergence of the stress-energy tensor to zero. This is the 
so-called "adiabatic" condition. We will see that as a result of this condition, one must take either G 
or A as an additional dynamical variable, and thus the dynamical system remains underdetermined. We 
would like to point out that if instead of retaining Einstein's Field Equations one retains the gravitational 
action, a new modified set of Varluminopic Field Equations will result, as we will see in Part IIIII Using 
the modified Field Equations, the adiabatic condition will completely determine the resulting dynamics 
without the necessity of additional dynamical variables. The considerations of this part of the paper are 
nonetheless interesting; and as the dynamical solutions in both parts are similar, the work done in this 
part will provide insight into the dynamics of the latter. 



6 Einstein's Field Equations 

6.1 The Metric Ansatz 

In this section we state the Christoffel Symbols and nonzero components of the Einstein and Ricci Tensors 
for the metric ansatz stated in 0. 

ChristofTel Symbols 

The nonzero Christoffel Symbols are: 

(36) 
(37) 
(38) 
(39) 
(40) 
(41) 

(42) 
(43) 
(44) 
(45) 
(46) 
(47) 



PP 


— a 


^ ee 


= a sin^ p 


"pa 

J- 00 


= a sin^ p sin^ 6 


ap 




pp 


= — sin p cos p 


pp 
^00 


~ — sin p sin^ 9 cos p 


^ a0 


= a-i 


i pe 


= cot p 


i 00 


= — sin 6 cos 6 


p0 

^ a0 


= a-' 


p0 

^ P0 


~ cot p 


p0 


= cote 
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The Ricci Tensor 

The nonzero components of the Ricci Tensor are: 



Rpp 
Ree 



^ 4 



4 sin^ p 

4 sin^ 9 sin^ p 



(48) 
(49) 
(50) 



The Ricci Scalar 

The Ricci Scalar becomes: 



The Einstein Tensor 

The nonzero components of the Einstein Tensor are: 



Gpp 

Gee 



_6_ 

-2 

—2 sin^ p 

—2 sin^ psin^ i 



(51) 



(52) 

(53) 
(54) 
(55) 



The Stress-Energy Tensor 

As is standard, we wiU be taking the stress-energy tensor to be that of a perfect fluid 

T^. = [e + P)U^U,- Pgp, 
where e is the energy density, P the pressure, and 

= ^ HI, 0,0,0) 



(56) 



(57) 



the four-velocity of a particle in the cosmological reference frame. The stress-energy tensor can be 
expressed in matrix form as follows: 



6.2 Einstein's Field Equations 











-Pg 








V 






\ 





-P.g00 / 



(58) 



Einstein's Field Equations (with cosmological constant A) can be written in cither of the following two 
forms: 



G p,ii — R^xu rfRg^-i^ 



;_4 J- fiu i^Sfiv 



R 



a 
SnG 



T^u - 7;Tgp.u - Ag, 



Einstein's Field Equations in the form of (|59|l imply the conditions 



2 



8ttG 
-8ttG 



P + A 



(59) 
(60) 

(61) 
(62) 
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Alternatively, the Field Equations written as l|6U|l imply the equivalent conditions 



For independent equations of motion, we will choose the following: 



= ^(s + 3P)-A (63) 
i = ^(e-P)+A (64) 







-A 


(- 


SttG 




d4 


(- 






SttG 





(65) 

In SjHl we will solve these equations explicitly for the case of a matter dominated universe, i.e. assuming 
E = eV is constant. 



We can define the following density parameters: 



Aa2 



= -J-^ (67) 



r!A = — (68) 
17fe = -fc (69) 

where fc = +1 in our model, as our metric Ansatz is equivalent to an fc = +1 FRW metric with the 
additional kinematic constraint c{t)dt = da. By convention, fig is the density parameter of the total 
energy density: matter, radiation, and any additional vacuum energy besides the cosmological constant. 
With these definitions, (|65|l becomes: 

ntot = n, + nA + nk = i (70) 

which is just the regular Friedman's Equation. We note that these definitions are comenserate with the 
standard ones if you simplify using the kinematic condition a = c. 

We would also like to note that Qg, = Qf/ + + Qy is taken to have contributions from radiation, 
matter, and nonlambdic vacuum energy; obtained by using £fl, sm, and ey into (|67(l . respectively. 

6.3 The Equation of State 

Since Einstein's Field Equations determine the energy density e and pressure P as a function of the scale 
factor, via (|65|l and (jHEJ, they therefore prescribe as well the equation of state. Introducing the state 
variable x = Aa^, we have: 

w{x) = - = (71) 



In particular, we observe 



6-x 



w{Q) = (72) 

w{2) = (73) 

^(3) = I (74) 

lim w{x) = +00 (75) 

lim u'(x) = — oo (76) 

lim w{x) = —1 (77) 

X^oo 



In particular, we note the following important epochs. 
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Equation of State 

6 I 1 1 1 



5 - 



4 - 



2 - 








12 3 4 5 6 

state variable: x = A a^ 



Figure 1: Equation of State vs. state variable x = 

6.3.1 Important Epochs 
The Null-x Epoch 

When X = 0, we have the equation of state P — —e/3. 
The Dust Epoch 

The universe can be treated as dust for x = 2- For x < 2, the pressure in the universe is negative. For 
X > 2, the pressure is positive. We are currently very close to this important epoch. 

The Radiation Epoch 

The universe can be treated as radiation for x = 3, when its equation of state is P = e/3. 
The Vanishing Epoch 

As X ^ 6^, the energy density goes to zero. At x = 6, the universe is completely empty in net energy 
content. After this epoch passes, the universe (if it continues) shall have negative energy density for all 
subsequent time. 

The x-Infinitum Epoch 

As X ~* ^1 the universe becomes dominated with a vacuum energy (unrelated to the cosmological 
constant term) with an ultimate equation of state P = —e. 

6.3.2 X and the the Cosmological Constant 

As we are considering a class of cosmologies where the speed of light is allowed to vary with cosmological 
time, there is nothing which prohibits the cosmological constant A from doing the same. For isolambdic 
cosmologies (with A = 0), x is an increasing function of cosmological time t, since a{t) is. This is not 
necessarily true for their varlambdic counterparts, where evolution of x is determined by x = 2Aaa+Aa^. 
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7 The Classical Thermodynamic Analogy 



7.1 Covariant Divergence of the Field Equations 

The Field Equations can be written in the form 



where 





£ 











\ 







-p 
















-p 







v 











-p 


) 



(78) 



(79) 



is the standard 1-1 stress-energy tensor. 

For the moment, we wiU aUow both G and A to vary with cosmological time, for the sake of generahty. 
We only include the case G ^ momentarily, to leave room for future speculations. Recognizing 



d_ 

da 



dt d 
da dt 



we obtain for the only nontrivially nonzero component of 

, 3(e 



V^(87rGT^^/d'^ + KSi^) the following 



SttG 



P) 



\ 87re 


(g 


AGd\ 


A 


)+- 




d5 ) 


+ - 
a 



(80) 



(81) 



This can be compactified into the following equation 

^ ^ iwa 



d 
dt 



In 



ea^G 



Thus 



3wd 



sa^G 



oc exp 



87rG£ 



87rGe 



dt 



(82) 



(83) 



We would like to point out that H81() imposes no additional constraint on the evolutionary dynamics of 
the system. By e and P satisfying H65|l and (|66|l . as solutions to the Field Equations, they autmoatically 
satisfy this condition. Whereas the standard FRW metric produces the Frieman Equation and the 
acceleration equation, and leaves open the equation of state; our modification, when submitted to the 
Field Equations, determines the equation of state and leaves open way for an additional constraint which 
one must impose on the system to determine the dynamics. For example, one now has room to impose 
the condition of isoergiticity on the universe, or adiabaticity, for example, as we will see in the next 
section. 



7.2 The First Law of Thermodynamics 

In this section we shall proceed entirely by analogy. Viewing the universe as a closed and isolated system, 
we can compare various terms which arise in (|81|) to terms which would arise doing a purely classical 
first law control volume analysis. 

Allowing E ~ eV be the total energy in the universe, where V = 27r'^a^ is the total volume of S^, 
we can apply the classical First Law of Thermodynamics and compare with our conservation condition 
dHII): 

E^Q-W^Q~pV (84) 
For the given E = 27r^a^e and V ~ 2'K^a^ , we have 

GTr^a^de + 27r^a^£ = Q - Qp-n'^a^a (85) 
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We can meanwhile rewrite (|81|l as follows: 



c 2 2- , o 2 3- / Sn^a^ed Ad'^a^n . iG] 22. /nc^ 

6n a ae + zn a e = — ztt ea — — 6p7r a a (86) 

\ a 4G G / 



Thus we can identify 



G 

— 27r'^ea'*— Cosmological Heating (87) 



4G G 



as a Cosmological Heating term, in analogy to classical Thermodynamics. 

In isoluminopic models, we have that (5 = 0. And thus it is a direct result of the form of Einstein's 
Field Equations that we must take V^T/f = (which account for the remaining terms of H86|) 'l. In a 
varluminopic model, if one uses the classical Field Equations, we are left with the necessity of imposing 
an additional constraint, which should take the form of a conservation law. We now mention two 
philosophicaly different approaches. 

7.2.1 The "Isoergic" Condition 

In a so-called isoergic model, we would take E = Q. The varying speed of light gives us room so that 
varluminopic effects (encapsulated in the Q term) can do the expansion work for us, so that the total 
physical energy of the universe is conserved. In this choice, physical energy cannot be converted into the 
gravitational energy of the universe, as is normally done in the isoluminopic case (though we note that 
in the isoluminopic case, this choice does not exist, and one is forced to consider a transfer of physical 
energy to a gravitational cosmological energy). The condition that i? = 0, however, is enough to provide 
us with an acceleration equation, so that the dynamics of the evolving universe can be analyzed. We 
will do this in §98. 

7.2.2 The "Adiabatic" Condition 

On the other hand, we do not have to abandon adiabaticity, which would impose the following dynamical 
equation of motion 

= |-27r^ea''— Adiabatic Condition (88) 

a 4G G 

This is tantamount to imposing the condition that V^T^ = 0. However, it is clear that if one is to obtain 
nontrivial dynamical solutions a ^ 0, one must include either G or A or both as dynamical variables. 
Thus the system remains underdetermined, and an additional constraint is needed to determine the 
dynamics. Such a constraint can be obtained by fixing the equation of state. We will explore this 
approach in iilHI 

It is interesting to note that this is not the case with the modified Field Equations we will derive 
in Part IIIII With the Varluminopic Field Equations one can impose the adiabatic condition V^T^ = 
and completely determine the resulting dynamics, without the necessity of the introduction of additional 
dynamic variables. 



8 Isoergic Models 

As mentioned previously, in these new models, we now have the freedom to constrain the total energy 
E ~ eV of the universe to be constant. In the isoergic model, the cosmological heating H87|l is what 
does the expansion work, so that the net physical energy in the universe is conserved. We will take 
G = A = 0. Using these conditions in H82|l leaves us with the following equations of motion 

ad 3w 

— = -r 89 



Recalling (|71|l and defining c = ct, we can rewrite this as 

dc _ -1 6- 3Aa2 
c ~ 4 6a - Aa3 



da (90) 
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Integrating, we find 

a = k{<6a- Ka^)-^/^ (91) 
Alternatively, we can proceed directly from the condition E = const.: 

STrGea^ = (6a - Aa^) = k'^ (92) 



where k — \J AGEtoi/'K. Solving for a wc arrive at: 

da k 



dt (6a-Aa3)i/4 



(93) 



which agrees with H91|l . but does not require that G and A be constant. Assuming G and A are constants, 
we can integrate H93(l to obtain 

4 if-^ ^ /-I 5 13 Aa2\ 



6 

where 2F1 is the Gauss Hypergeometric Function (see |18j): 



.F,(a,/?;X<5) = EW§''^" (^5) 

n=0 



which has radius of convergence \5\ < 1. Here, (a)„ is the shifted factorial 



, , T(a + n) , , 

(«)« = (96) 

and T{z) is the T function: 

T{z) = / t^-ig-* dt (97) 
Jo 

The solution can be simplified as: 



kt : 



4^6^ >^r(n-l/4) /Aa2^" 



1/4) ^„ (5 + 8n)nl \ 6 / ^ ^ 



r(-l/4) ^ (5 + 8n)ri! V 6 
We plot the speed of light in Fig. [3 left and the scale factor in Fig. 13 right. 




Figure 2: a(i)/a(to) vs. i (left) and c(t)/c(to) vs. t (right) 

The radius of convergence of H98|) is Aa^ < 6. At the epoch Aa^ = 6, the total energy density of the 
universe vanishes. Beyond this point, the energy density would have to become negative for the universe 
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to continue, a situation we view as unphysical. We can use H94f) with the following mathematical fact 
(which holds for c — — a > 0, as in our case): 



iFi{a,b] c; 1) 



T{c)T{c-b-a) 
T{c-b)T{c-a) 



to see that this model predicts an end to the universe in the finite time 

_ 4.6V«r(13/8)r(5/4) ^ 



(99) 



(100) 



where we used a = 1.32 x 10^6m and Aa^ = 2.1 (so that JIa = .7). These numbers imply the present 
age of the universe to be close to 12 Gyr, so that we are currently 8.4 Gyr from the end of the universe. 

Moreover, using the redshift formula H25|l . the luminosity distance formula (|34|l . and the definition 
(|35|l . we can compare the theoretical curve produced by this model against data points from Supernova 
la data as recorded by Riess ^S]- We find our theoretical curve matches the supernova quite well, as is 
shown in Fig. O 



Supernovas la Data. Riess 1 998, x: Nearby SN, o: High-z SN, — : theory 



^40 - 
^J8 - 
36 - 
34 - 




Figure 3: Supernovac la data vs. theoretical predictions 



Without the Comsological Constant A, this solution degenerates to 

1 /5fct\"/' 



(101) 



Inverting gives the scale factor, a{t), as a function of time. 



Part III 

Varluminopic Field Theory 



In Isoluminopic Field Theory, Einstein's Field Equations can be derived by varying the following action 

16^G 



5"= / {-^{D\ + 2Ail)-Tl\dtd^e (102) 



Gothicized variables denote tensor densities. For instance, = R^—g is the tensor density of the 
Ricci scalar, and it = ^/— 5 is the tensor density of unity. 9Jl represents the Lagrangian density of the 
mass-energy which occupies the space-time, eg. matter, radiation, etc. It is well known that one cannot 
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construct a scalar density form the metric and its first derivatives alone. Thus the Ricci scalar density 
represents the most general covariant candidate for a gravitational action. Even though it contains 
second derivatives of the metric, these have no effect when varying the action ()102|l . see Landau and 
Lifshitz 13, Tolman or Pauli ^I], amongst others. Following the derivation in Landau and Lifshitz 
with appropriate modifications, we will show in ij^lthat this is true even when taking the time- varying 
speed of light into account. 

The action for any Varluminopic Field Theory should therefore reduce to ()102ll in the isoluminopic 
limit. The most natural candidate for an action is therefore (|l()2(l itself. In HIOI we will derive the 
resulting field equations from the action p()2|l . taking the varying speed of light into account. We will 
then apply these new field equations to the k ~ +1 FRW model, including the kinematic constraint d = c, 
to determine the dynamical equations of motion of the resulting theory. We will see that imposition of 
the additional energy conservation constraint V^T^ = is enough to determine the resulting dynamics 
explicitly; unlike with the classical Field Equations, where one was left with the necessity of having to 
also include G or A as dynamical variables. 

In formulating the variational principle for gravity, we will find that the cosmological coordinates 
{t, p, 9, 4>) will be better suited for our needs. We will keep these coordinates in the back of our mind 
while we derive the modified field equations. There are two important features in doing this. Notice 
that all factors of c in H102() are attached solely to the gravitational action. Thus, by letting a;° = t, as 
is in our prerogative to do so, we decouple the varluminopic effects from the matter action. In standard 
classical derivations of the Field Equations, the entire integrand of Hl()2|) is typically divided by an extra 
factor of c, in exchange for swapping the dt with a da (in the notation allowed by our insight da = cdt). 
Moreover, in this choice of variables we have 900 = c^, so that variations of c can be easily related to 
variations of the metric, as will be shown in HIOI 



9 The Action of the Gravitational Field 

As is standard in relativity theory, one identifies the space-time metric with the graviational field. It 
is well-known that there is no scalar density which depends only on the metric and its first derivatives. 
Thus the Ricci tensor is taken as the key ingredient of the Lagrangian density of the gravitational field. 
Including factors of c which must be present for our later considerations, the action is similar to the 
nominal action 



^mtd-^e 



where 



is the Ricci scalar density. Our aim is to show that this covariant Lagrangian density is equivalent to 
one which only involves the metric and its first derivatives (but depends on the coordinates). Thus we 
will show that 

dm' 



j c^mtd^e^ J (tdtd^e + J ^dtd^o 



where to* is a vector density, which, by means of Gauss' Theorem, can be converted to an integral 
over the boundary and therefore ignored, and 2; is the pseudo-scalar density which depends only on the 
metric and its first derivative. We will follow closely the discussion in Landau and Lifschitz [7] §93, with 
appropriate modification for the varluminopic effects. 
The integrand may be expanded as 



-9\9 



dx^ 



9 i ik^ Im 



ikT^mj^l 
9 i ^ km 



(103) 



Notice the first two terms may be written as 



-99 



-99 



■ik 



dAk 


d 


dx^ 


dx'- 


Ml 


d 


dx'^ 


dx'^ 



dx'' 



"99 



ik\ 



(104) 
(105) 
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The first term on the right hand side of either of these equations may be dropped as it may be converted 
into a boundary integral which vanishes, as we wiU take the variations to vanish on the boundary. Thus 
we have 



+^ZV^99''"-^ - (106) 

4: ik fY<7nj^l j^l j^m \ , j^m „ik ) j^l ^jk ) 

— il ^ km ~ ^ ik^ Im) ^ iniQ Q^l ' ) 

This pseudo-scalar density only depends on the metric and its first derivative, and we have that 

6 J c^mtd^e = 5 j (Bdtd^e (i09) 

Thus we shall not hesitate in using [H as a Lagrangian density, even though it contains second derivatives 
of the metric, as they have no effect when one takes the variation. 

10 The Varluminopic Field Equations 

We wish to find the corresponding equations of motion by varying the action 1)102(1 with respect to the 
space-time metric. The only difference between our procedure and the classical approach is that c now 
varies with cosmological time. Hence we need to determine the effects of this on the resulting equations 
of motion. Noting that 

SV^J = -lV^9^kSg''' (110) 
we see that the variation of the first part of the integrand of H102|) is 

S (c'^m) = S {c^g'''R,kV^ = (^^,k - ^'y^g^k^ Sg''' + c^Q'^SR.k + 49^c35c (111) 

In order to deal with the c^Q^^SRik term, we will again closely follow the derivation in Landau and 
Liftshitz with the appropriate modifications necessary for our varluminopic considerations. Choosing 
a locally geodesic frame, we have that 

d'g^HR.k - c^g'^-^(<5rL) - c4/^(<5r^,) (112) 

Now consider the vector 

= c^g^Hv\^ - c^g^'STt^ (113) 
Taking the divergence with respect to a;', we find 

9?"^ d^^ 9^('^r,J-H(5 ST.^k-g ST^k)—^- (114) 

Returning to an arbitrary reference frame, we see that 

^9''SR.^ = 4=S - C^'^rL - /fff,) ^ (115) 



-Q^l V.V -^k y g^l 

We will now show that the extra term which arises (second term on the right hand side) actually vanishes. 
Recalling that 

Sgip = -~gikgps5g^' (116) 
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we have that 



6T: 



ik 



5 (5/./-^r,L) 



■ik 



91 



Thus 

and therefore 
so that 



ik 



g^'^6T\,-g^^ST^,^0 



dtcPe 



an 



(117) 
(118) 
(119) 
(120) 

(121) 
(122) 

(123) 



as we take all variations to vanish on the boundary. 

We are therefore justified to write the variation of the full action (|102|l 



5 



167rG 



(SH + 2 All) - Wl \ dtd^e 



($n + 2AiI) Sc 



dtd^e 



where 



25,: 



1 



ik 



(124) 



(125) 



is the Einstein tensor density. 

Writing the variation of (|1()2|1 in the form H124() , we are now able to see the full advantage of consid- 
ering our variational principle using the coordinate x*^ = t as cosmological time as opposed to = a, as 
is classically the choice. First, one no longer has a factor of appearing in the matter action, so that 
the varluminopic effects limit themselves to effects on the action of the gravitational field. Moreover, 
with the choice of t for our cosmological coordinates, we have the relation 



9 



00 



so that 



or 



5g 



00 



7 5c 



2 ^ 



(126) 

(127) 
(128) 



and thus the variation of 11U2|I can be written, for our choice of cosmological coordinates, as: 



167rG 



(fR + 2AU) - 9Jt \ dtd^e 



^(6,,-A0,,)-^ 



SnG 



(m + 2Ail) S'^S'i } Sg^'dtd 



Thus, the modified Field Equations, in our choice of cosmological coordinates, become 



(129) 



(130) 



Consider now a general coordinate system x^^. To every point in space-time we identify a cosmological 
scale factor a(x'^), which represents a scalar function on the space-time. As the map a{t) is one-to-one, 
which identifies the cosmological age of the universe with the scale factor, we also have the scalar function 
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which identifies the cosmological age of the universe to every point in the space time. Transforming 
(|130|1 to an arbitrary choice of coordinates, we therefore find 

G,. = ^T,. + Kg,^ + 2c\R + 2A)^^ (131) 

Moreover, recalhng da = c{t)dt, we can ahernatively express the extra term as a tensor product of the 
divergence of the scale factor a{x^) with itself: 

^ SttG ^ , da da ,,„„^ 

G,. = —T,^ + A.g,. + 2{R + 2A) — — (132) 

Hence, even though we choose to work in the cosmological frame to simplify computations, the modified 
Field Equations can nonetheless be written in an entirely covariant fashion. 

11 A Cosmological Action 

Suppose we would like to add a cosmological action of the form 

Sc= [ (tdtd^e (133) 



to the action p()2(l , where £ is of the form 

€^C{a,c)il (134) 

We have 

dC dC 

d(f = —USa + -T^iiSc + CSii (135) 
da ac 

However, noting that 

Sc^'^^Sa (136) 
dtda ^ ' 

we can write 5a in terms of a corresponding 5c. Using this relation with l|ll()(l and 1)128(1 . we find 

«-i4f?^ + ^VX-%^|vS;*.- (137) 



2 \c da dc 

Adding this to ((129() and converting to an arbitrary coordinate system, we obtain the following field 
equations: 

SttG /. C\ / . ^TiG fade dC\\ da da , , 

However, for the rest of the discussion here we will take £ = 0, so that our action is to coincide with the 
classical action for general relativity. 

12 Application to our Varluminopic A; = +1 Model 

Using the Varluminopic Field Equations H132|) with our metric Ansatz © and the stress-energy tensor for 
a perfect fiuid 15()|) . coupled with our hypothesis 10} , we obtain the following Field Equations (compare 
with ^ and (EH): 

6 SttG , 24 

= ^^£ + A- — +4A 139) 
a^ a^ a^ 

4 = -^^ + A ,140, 
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These may alternatively be written in the form (compare with (|65|l and l|6tif) ): 



5(2* 

8^ 

i4 



P = 



a' / . 2 



8ttG 



(141) 
(142) 



These have the same qualitative features as those for the classical Field Equations approach as discussed 
In particular, notice that (|141|l implies a new Friedmann Equation: 



fie + S^A + Silfc = 1 



Also, the equation of state is now given by 



(143) 



(144) 



where x — (compare with 1)71(1 V The equation of state is plotted in Fig. 0] We also included the 
equation of state obtained by using the classical field equations ((71|l on the same plot for reference. The 
equation of state for the modified field equations, as given by 1)144(1 . is the curve which has been vertically 
compressed by a factor of 5. 



Equation of State 




state variable: x = A a 



Figure 4: Equation of State vs. state variable x — ^o? for classical and modified Field Equations 

In particular, we note that the matter-dominated epoch occurs at x = Ka? = 2 and that the radiation- 
dominated epoch occurs at x = Aa^ = 4.5. Recall that using the classical field equations, as in these 
epochs occured at x = 2 and x = 3, respectively. 

Applying the additional constraint that ^ ^iTj^ = 0, one can derive (either directly from the stress- 
energy tensor, or by incorporating this condition when taking the covariant divergence of the Varlu- 
mmopic Field Equations ^^): 

ad 



Setting c = a, we have 

1 dc 
c da 

And therefore, integrating, we obtain: 

da 
It ~ 



1 6 - 3Aa2 
5 6- Aa2 

1 6 - 3Aa2 
"56a- Aa3 



(6a- Aa3)i/5 



(145) 
(146) 



(147) 
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where fc is a constant of integration, different than the constant of integration which appears in (|93|1 . 
(|147|l may now be integrated to find 



This can be simphfied to 



H^^^6.S^(4,||^) (148) 



2r{^l/5)^ (3 + 5n)n\ [ 6 ) ^'^'^ 



The scale factor ()149|l and the speed of hght (|147|l are plotted against time in Fig. We would like 
to point out that, despite appearences, we are slightly past the minimum value Cmm in Fig. ()147|l r. We 
used a(to) = 1.32 x lO^^m and Aa^ = 2.1 (so that f^A = -7). 



Evolution of the Scale Factor 



Evolution of the Speed of Light 




Age of the Universe Gyr 



Age of the Universe Gyr 



Figure 5: a{t)/a{to) vs. t (left) and c(^)/c(^o) vs. t (right) 

Moreover, using the redshift formula H25fl . the luminosity distance formula (|34fl . and the definition 
(j35|l . we can compare the theoretical curve produced by this model against data points from Supernova 
la data as recorded by Riess jl3| . We find our theoretical curve matches the supernova quite well, as is 
shown in Fig. El 

Supernovae la Data, Riess 1 998, x: Nearby SN, o: High-z SN, — : theory 

48 
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44 
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36 
34 

32 
30 
28 

10"^ 10"^ 10"^ lo" 10^ 

redshift z 




Figure 6: Supernova la data vs. theoretical prediction 
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We would like to note that the scale factor evolution, as seen in Fig. O and the theoretical redshift- 
distance curve, Fig. (BJ are nearly identical in appearance to those predicted by the isoergic model from 
3H1 See Fig. Eland Fig. El 

At the epoch Ka? = 6, we see that the speed of light approaches infinity. This epoch also corresponds 
to the dissapearance of mass-energy, as is indicated by the density equation H141(l . We interpret this as 
an indication that the universe has "used up" the total amount of mass-energy in the universe by doing 
the expansion work. This is essentially the death of the universe. If one were to continue Ka? beyond 
the value 6, one would obtain negative energy densities, which we view to be unphysical. Referring to 
(j99|l . we can determine the cosmological time of the universe's impending demise to be: 

This model predicts that the current age of the universe is approximately 12.3 Gyr, so that the universe 
would have about 8.7 Gyr of life left. 

Alternatively, without a cosmological constant, the solution to (|147|l would reduce to 

13 A Note About the Varlambdic Approach 

By introducing a time- varying cosmological constant, one can effectively "control" the equations of state 
as given by the Classical Field Theory (CFT) and Varluminopic Field Theory (VFT) approaches. The 
radiation- and matter-dominated epochs occur at 

XZI=XIII = 2 (152) 
X?ad^ = 3 (153) 
xYId" = 4.5 (154) 

where x is the state variable x — Aa^, see (|71|l for the CFT case and (jl44|l for the VFT case. Replacing 
A with X as a dynamical variable, we can choose a function x{t) which is initially constant with Xrad 
and whose value changes to Xmat during a rapid phase transition. This corresponds to a rapid phase 
transition from a radiation- to matter-dominated universe, with time varying cosmological constant. 
The Field Equations can be rewritten as 

see (ESI, (ESI (CFT case) and ifTITll . (tH^ (VFT case). For both cases, the conservation law V^T^ = 
can be written as 

|+3(e + P)a = (158) 



a 



Defining the acceleration parameter 
we have that p58(l produces 



g = ^ (159) 



„ 1 ay — 12d 

1 ax-12«(7-x)/5 
^^"^ - 2+ M6^) 



22 



Interestingly enough, recalling l|152|) - (|154|l . we see that both approaches lead to the same steady-equation- 
of-state dynamics: 

Qmat = -J (162) 
Qrad = (163) 

Thus, for either the CFT or VFT approach, the scale factor becomes 

a(i)rad (X t'/^ (165) 
for a matter-dominated and radiation-dominated universe, respectively. 



Part IV 

Conclusion 

14 Conclusion 

In this paper, we presented a new cosmology where the speed of light varies with cosmological time 
subject to a fundamental constraint (0J. We explored the implications of this in Part |n] using the 
classical field equations and showed that it leads to a new and interesting dynamics of the scale factor. 
In Part IIIII we took a different approach. Retaining the classical gravitational action (|102(l . we showed 
that one obtains a new set of Field Equations (|132(l . With the additional constraint V^T/;* = 0, we solved 
these modified field equations explicitly p49|l and found that they imply a possible end to the universe 
which would occur at the epoch Aa^ = 6, or roughly 9 Gyr from now. The model predicts a universe 
which is currently around 12 Gyr old, which does not contradict known timelines. A key feature of the 
model is that it solves the Horizon Problem without the need of inflation. 

We showed that, although the redshift formula (|25|l remains unaltered, the formula for the luminosity 
distance H34|l inherits an extra factor in the varluminopic case. Using the dynamical solutions as presented 
in fJHl and iil2l we plot the redshift-distance curve and showed that in both cases they match with the 
experimental measurements recorded in Riess ^3], see Fig. O and Fig. |H| These results need to be 
examined more closely, but the preliminary figures included in this paper are a good sign. 

Big Bang nucleosynthesis, on the other hand, could lead to problems. Our model changes the con- 
ditions at the time of nucleosynthesis quite drastically. When the universe was at the temperature T = 
1 MeV, the speed of light, given by H147|l . was 229 times larger than its present day value. Thus it is 
not immediately obvious whether nucleosynthesis will still work in our model, and will considered in 
future works. Other astrophysical issues, such as the cosmic microwave background and the formation 
of structure, also need to be studied if this model is to remain viable. 

Our model changes the history of the early universe. We no longer have the need for inflation. Phase 
transitions should work out to be different, eg. electroweak and quark hadron. And further it changes 
the conditions during baryogenesis, which might also be different. All of these issues are uncertain in the 
standard model, and they may be better or worse here. This will be the focus of much further research. 
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